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ABSTRACT 


Predicting the effectiveness of an aerially sprayed pesticide 
requires a method for calculating pesticide concentration and the 
distribution and temporal variation of the size of the droplets 
comprising the spray. The rate of evaporation of the spray 
droplets is a critical factor affecting these calculations. This 
report describes some methods for calculating the evaporation 
rate of a droplet, with particular emphasis on how the governing 
equations were derived. These calculations will enable those 
working in aerial application to compare data from other experi- 
ments without having to reproduce the tests to evaluate the 
results. 


A report on Equipment Development and Test Project 2664, Optimum 
Size Spray Droplets, funded by the Forest Insect and Disease 
Management Staff. 


PESTICIDE PRECAUTIONARY STATEMENT 


This publication reports research involving pesticides. It does not contain 
recommendations for their use, nor does it imply that the uses discussed here 
have been registered. All uses of pesticides must be registered by appropriate 
State and/or Federal agencies before they can be recommended. 


CAUTION: Pesticides can be injurious to humans, domestic animals, 
desirable plants, and fish or other wildlife — if they are not handled or 
applied properly. Use all pesticides selectively and carefully. Follow 


recommended practices for the disposal of surplus pesticides and pesticide 
containers. 
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INTRODUCTION 


An important objective of the Forest 
Service aerial applications program is 
the development of a technique for pre- 
dicting the behavior of a pesticide 
spray once it has been released from 

the aircraft. This prediction requires 
a description of the temporal and 
spatial distribution of the droplets 
comprising the spray. A critical factor 
affecting this prediction is the evapor- 
ation rate of the spray droplets. 
Because droplet evaporation data are 
exceedingly difficult to obtain, the 
Missoula Equipment Development Center 
(MEDC) has contracted with an outside 
agency to collect data from evaporation 
experiments. To assess the validity of 
these data, Center engineers investigated 
several theoretical models, summarized 
in this report, for predicting droplet 
evaporation rates. 


No satisfactory theoretical model yet 
exists for predicting the evaporation 
rates of multiphase liquid droplets or 
droplets containing particulate matter. 
Experiments with these types of droplets 
have yielded only scattered results. 
Consequently, calculations for evapora- 
tion rates in this report have been 
based on pure liquid droplets so that 
the results can serve as a standard 

for assessing evaporation rate data 
obtained from actual experiments. 


Two theoretical models for predicting 
the evaporation rates of pure liquid 
droplets have been found to produce 
results that are in close agreement with 
results obtained from droplet evapora- 
tion experiments. These models, based 
on the works of J.C. Maxwell, 

I. Langmuir, and N.A. Fuchs, are summar- 
ized in this report. For each model, 
the assumptions used to derive the 
governing equations have been carefully 
stated and the range of applicability 
specified. 


In addition to excluding multiphase 
liquid droplets and particulate droplets 
from the theoretical calculations for 
droplet evaporation rates, the effects 
of condensation and the influence of 
gaseous droplets adjacent to the spray 
cloud also have been excluded. 


Evaporation is the process where mole- 
cules of a liquid leave the liquid and 
join the surrounding gaseous or vacuous 
state. It also may be the case that 
some of the vapor molecules of the 
surrounding medium collide with the drop- 
let surface and thus join the liquid. 
This process is termed condensation. In 
this report, it will be assumed that the 
process of condensation may be ignored, 
however, this assumption may not be 
valid for all aerial spray applications. 
For example, if the spray carrier were 
water-based and the relative humidity 

of the air were high, a thorough 
discussion of droplet evaporation would 
have to include the effects of condensa- 
tion. Condensation also may be important 
in spray clouds that are very dense 
because evaporation from one set of 
Cloud particles will contribute to con- 
densation on another set of cloud 
droplets. Nevertheless, in this report, 
emphasis will be given to examining 
expressions for the evaporation rate of 
single droplets and condensation 

effects will not be considered. 


From the laws of hydrodynamics and 
thermodynamics it is known that pressure 
and temperature gradients in both the 
surrounding gas and the droplet govern 
the motion of the molecules comprising 
the gas and the liquid droplets. A 
discussion of the mass transfer processes 
induced by these gradients must be 
derived from a statistical mechanical 
description of the behavior of matter. 
This description has not yet been com- 
pletely worked out. Consequently, this 
report will adopt the macroscopic or 
continuum mechanics approach rather than 
the particle or statistical mechanics 
approach. 


Droplet evaporation is only one of many 
factors that affects the behavior of 
pesticide sprays. Other phenomena, such 
as topography, wind velocity, and turbu- 
lence are equally important. Thus, the 
theoretical models in this report provide 
a satisfactory "first approximation" for 
accounting for the affect of droplet 
evaporation on the behavior of aerial 
sprays. 


The first portion of this work will 
present the development of expressions 


for determining the evaporation rates 

of droplets that are at rest with respect 
to the surrounding medium. In the 

second portion of the report, expressions 
for the evaporation rates of droplets 
that are moving with respect to the 
surrounding medium will then be obtained 
by modifying the previously developed 

"at rest" results. 


An excellent survey and critique of 
experimental methods for determining 
evaporation rates has been given by 

Beard and Pruppacher (1971). Their paper 
also presents useful equations for pure 
water droplets. Pruppacher and Rasmussen 
(1979) have recently extended the pre- 
vious work to include the determination 
of the rate of evaporation of large water 
drops falling at terminal velocity in 
air. Woo and Hamielec (1971) present an 
accurate numerical method determining 
evaporation rates and discuss the work 

of other researchers. 


EVAPORATION OF STATIONARY DROPLETS 


Maxwell-Langmuir Model of Evaporation 


The earliest and simplest model of evap- 
oration is Maxwell's (1890). His results 
were used by Langmuir (1918) to obtain an 
equation governing the time rate of 
change of the mass of an evaporating 
droplet at rest with respect to the 
surrounding medium. The Langmuir 
equation may be used to derive an equa- 
tion for the temporal variation of the 
droplet mass, surface area, and radius. 
To enable the reader to more easily 
follow the derivation of the model 
equations, we first give a brief outline 
of the development. 


1. The evaporating droplet is 
assumed to produce an infinitely thin 
shell of vapor molecules which then move 
or diffuse out into the surrounding 
medium. The concentration of molecules 
in the shell is assumed to be at the 
saturation point. 


2. From assumptions governing the 
medium and the shape of the droplet 
during the evaporation process, an 
expression is derived for the rate at 
which mass is transferred from the shell 
to the surrounding medium. 


3. This expression is equated to 
the time rate of change of the droplet 
mass to give Langmuir's equation. 


4. lLangmuir's equation is then 
integrated to give the temporal variation 
of the droplet surface area, mass, and 
radius. 


Development of the Equations 
Assumptions: 


@ The 
remains so 
process. 

@® The droplet is motionless with 
respect to the infinite medium assumed 
to surround the droplet. 

e The process of evaporation results 
in a concentration or density of evapor- 
rating molecules in the spherical shell 
immediately surrounding the droplet and 
this concentration is equal to the vapor 
concentration of a saturated vapor 
corresponding to the temperature of the 
droplet. 

@ Diffusion is the mechanism 
governing the transport of the molecules 
in the spherical shell to the surround- 
ing medium, 


droplet is spherical and 
during the entire evaporation 


Derivation: 


The mixing of two substances by the 
passage of molecules from one substance 
through spaces between the molecules of 

a second substance is called diffusion. 
Diffusion is a mass transfer process 
Since a mass of molecules leaves one 
substance to intermingle or to be trans- 
ferred to the space occupied by the mole- 
cules of a second substance. This 
suggests that an important variable in 
describing a diffusion process is the 
concentration of molecules, that is, 

the number of molecules per unit volume. 
It will be postulated that molecules 

move from a region of higher concentra- 
tion to a region of lower concentration 
at a rate which is proportional to the 
difference in concentrations. Let C 
denote the con¢gentration of vapor mole- 
cules and let J denote the mass of vapor 
molecules passing through a unit area 
normal to the vapor mass flux ig uni 


ie 
: ; ; = 
time. J is expressed in gm cm Secu. 


Now, for a spherical surface of radius 
x (cm) 


> > 
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is the mass flux (gm exerts) leaving the 
surface through an element of surface 
area, ds, in a direction perpendicular 
to ds. The symbol ds denotes a.vector 
whose magnitude is the area (cm) of 
the small increment of the surface 
through which the particles are passing 
and whose direction is perpendicular 

to the small area increment. Because 
the surface is assumed to be spherical 
and to retain its spherical shape 
throughout the evaporation process, if 
spherical coordinates r, 9, and ¢ are 
introduced, the magnitude of an element 
of surface area is given by 


r2 sind dO dp ; 
The mass flux of vapor molecules through 
this area is then 
rae oe acy an 
Jmomnex “sind do’ ao 


Aa 


where n is the unit direction vector 
normal to the surface element, ds. 
Hence, the mass flux, F, of vapor 
through the entire spherical surface is 


6=1/2 p=2T ae 
F mS * n xr* sin dO@ do 6 


=-1/2 =o (1) 


Cay 


Because the vapor flux of molecules is 
assumed to be directly proportional to 
the gradient of the vapor concentration, 
Caaterollowss that 


ee Ware (2) 


where V is the gradient operator and D 
is the constant of proportionality or 
diffusion_coefficient. 2P is_jeasured 
in gm sec cm’ = gm_gm sec and C 

is measured in gm cm. It is worth 
noting that other investigators such as 
Green and Lane (1964), ‘Bradley, Evans 
and Whitlaw-Gray (1946) and Romig and 
Long (1969) express the vapor cgncen- 
tration, €, in units of mol cm rather 
than in units of gm cm” as is done in 
this work. Our’ Cis actually the 
density of the gas. 


Now, 


> “an a 
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and since, for a spherical surface, n 
is the unit outward normal vector in 
the direction of the radius, it follows 
that 


Hence, combining equations (1) and (2) 
yields 


= 71/2 27 


-Dr2 —= 


F = 5 


sin@ 46 dd 


CN 2 
Il 


aC 
ug 


oo 
ll 
ll 


-7™/2 6 fe) 


which, upon carrying out the integration, 
gives 
2 3C 
F = -47 Dr er (3) 
This is the equation used by Fuchs (1959) 
in beginning his work on droplet evapo- 
ration. 


It is now assumed that the process of 
evaporation is stationary, that is, the 
mass flux of vapor molecules across an 
element of the droplet surface is 
constant. This assumption is valid 
providing the concentration of vapor 


molecules in the far surrounding medium 


is much less than the concentration of 
molecules in the droplet. This is 
certainly true for droplets whose 
radius is larger than 1 micrometer. 

The assumption also requires there be 
no significant thermal gradients in the 


system. Thus, equation (3) can be 
written as 

Py ar = -4m D dC 

a 


which, upon integration, yields 


F -_ 
pp at DCs Cc 
where C is a constant of integration. 
For very large values of r it is assumed 
that the vapor concentration is the 
ambient vapor concentration in the 
surrounding medium, C,.. Thus, for 


Yr->o it is the case that C > C., 


whence, 
CimaATDC sy 5° 
Therefore, 
ly 
ms ins 4nD(C-C_) 
or 
Pom 4icD(C-C) (4) 


This is Maxwell's equation and it gives 
the total mass flux of vapor molecules 
leaving the imaginary shell of vapor 
surrounding the droplet to diffuse into 
the surrounding medium. The flux is 
given in terms of the diffusion coeffi- 
cient of the vapor and the difference 
in vapor concentration in the surface 
layer of the droplet and in the surround- 
ing medium. By assuming that the 
vapor obeys the ideal gas law 


Pm 
c= — : 
RT (5) 
equation (4) may be written as 
_ AnDmr 
ee ee ead (6) 


which is another form of Maxwell's 
equation. In equation (5), p is 
expressed in mm of Hg, T is expressed 

in degrees Kelvin and m is expressed in 
gm mol R is the universal gas constant 
whose value is 


R = 62,360 em 7aeg K/mol 


when the pressure is expressed in mm 
Hg. 


Let M denote the mass (gm) of the 
droplet. The equation governing the 
rate of evaporation of the droplet is 
then 


OM _ 5 

de = (7) 
Hence, 

dM 

— = D (C-C 

= AnxrD ( a) (8) 
or in terms of the vapor pressure, the 


equation is 
47Dmr 


Seen eee) 


RT (9) 


These are the two most common forms of 
Maxwell's and Langmuir's equations for 
the temporal change of eo Ee of 

the droplet. 


To obtain expressions for the temporal 
variation of the surface area, radius 
and mass of the droplet we begin by 
deriving the relevant differential 
equations. Since the droplet is assumed 
to be spherical, we have 


where p is the density (in gm on] of 
the droplet. Under the assumption that 
the droplet density is constant in 


time, it is the case that 
aM _ 9 ar 
dead ade 


The surface area, S, of a sphere is 


S = 4nr2 
and thus, 
ase sees fae 
at Ever 
- 2 aM 
oz.at (10) 
By using equation (10) with equations 
(8) and (9) we have 
OS is -87Dm cere) 
dt Smaps (11) 
or 
ds 2 -87Dm eee 
at prt ‘P Po . (12) 


By integration we obtain 


SED aecenyery C) 


‘and 


= 


(p-p,)t + c, 


To evaluate the constants of integration, 


Cy and C 2" we note that when t = ty = 0, 
S = §. = 4nr2 
fe) fe) 
where r_ is the radius of the droplet 
at the beginning of the evaporation 
process. Therefore, 


Cc. = 4nr? = c 
fe) 


and the equation governing the temporal 
variation of the surface area is 


“oary Sy WORD h4 
S 4nre = (C Co. EG (13) 
or 
Ss = 4nr2 - 87Dm 
ORT (p ts eam ae 


Since the droplet is assumed to be 
spherical during the entire evaporation 
process equations (13) and (14) may 
also be written as 


2 2 2D 
= pe = t 
E Shy (C-C,) (15) 
and 
=p 2 a Dm = 
See ort 'P p,,)t = (LO) 


These are the equations obtained by 
Maxwell and Langmuir for the temporal 
variation of the droplet radius. 


Discussion of Maxwell and Langmuir Equations 


Equations (13) and (14) show that the 
droplet surface area decreases linearly 
with time. Fuchs (1959) developed an 
alternative form of equation (14) to 
account for the fact that near an 
evaporating surface there must be a 
mass flow of vapor molecules directed 
away from the surface of the droplet. 
This flow is governed by the laws of 
hydrodynamics. Fuchs also showed that 
the correction to be applied to equation 
(14) to account for the mass flow is on 
the order of 1 percent or less for 
particles larger than 5 micrometers. 
For most insecticides, droplets whose 
diameters are less than 5 micrometers 
contain insufficient toxin to affect 
the target insect. Consequently, when 
including the effect of evaporation in 
estimating the probability of kill of a 
target insect population, it is not 
necessary to include Fuchs' modifica- 
tions of the equations of Langmuir and 
Maxwell. However, when estimating the 
environmental impact of the drift of 
aerial sprays, it will be necessary to 
account for the ultimate deposition of 
very small droplets. Such an accounting 
will require a modification of the 
evaporation equations of Langmuir and 
Maxwell. 


It is also possible to alter the previous 
development to account for the nonspheri- 
city of the evaporating droplet. However, 
it is thought that nonspherically 

shaped droplets result from an extremely 
complex interaction of such phenomena 

as high-speed shear flow, turbulence, 

the break up of the spray injection 
stream, variation in surface tension, 


etc. Since the analysis of this inter- 
action has not yet been worked out, the 
evaporation of nonspherical droplets 
will not be considered. Furthermore, 
Pruppacher and Rasmussen (1979) found 
that changes in droplet shape, internal 
circulation within the droplet, and 
droplet oscillation did not reflect 
themsleves in changes in the evaporation 
behavior of the droplet. 


It is reasonable to assume, that spray 
droplets sufficiently far removed from 
the aircraft spray nozzle become and 
remain spherical due to the surface 
tension of the liquid. Consequently, 
the analysis presented in this report 
should be applicable to such sprays. 


Equations (15) and (16) permit the 
calculation of the time, t_, necessary 
for complete evaporation of the droplet. 
Thus, setting r = 0, gives 


2 Zz 
IU Ure 
pre pRTr® 


houeop=(ceewy’ opmip=p.) | + <27) 


To illustrate the use of equation (17) 

we compare the evaporation times of 

water droplets of varying initial sizes 
under different climatological conditions. 
Consider the following cases: 


Dry bulb temperature Humidity 
Case surrounding medium (p/p x 100) 
(Gane) (percent) 
a 50 50 
b 50 60 
Sc 60 50 
da 60 60 


For water we use the following values 
for the constants appearing in equation 
(Ze 


6) 
p = 1 gm cm a 
D = 0.240 cm ec 
m = 18 gm mol, 
R = 62,360 cm /deg K/mol . 


The calculation is readily followed with the aid of table 1: 


Table 1.--Corresponding sets of values for equatton (17) 


Vapor 
Dry Wet Vapor Temp. pressure -5 
bulb bulb saturation surrounding surrounding peRTExedO 
Case temp. Humidity temp. pressure medium medium 2 Dm (p-p,) 
(" F) (percent) (— F) (mm Hg) CK) (mm Hg) 
a 50 50 41.9 6.769 283 3.384 6.036 
b 50 60 43.6 7.280 283 4.368 7.014 
c 60 50 5073 9.309 288.556 4.655 4.418 
d 60 60 52.4 10.058 288.556 6.035 Diehl 7 


Now, with the aid of equation (17), the following table for evaporation times 
is obtained: 


Table 2.--Evaporatton times of varying stze water droplets 
for different ratios of the surrounding temperature to 
the relattve humidity 


= (sec) 
ro (um) Case (a) Case (b) Case (c) Case (d) 
(50/50) (50/60) (60/50) (60/60) 
1000 6036 7014 yb Ilys 4418 
100 60.36 70.14 Ayia 44.18 
10 OSSI7 7 0.7014 OS5157:7 0.4418 


It is interesting to examine the consist- 
ency of these results. For example, 
suppose the relative humidity of the 
surrounding medium is held constant and 
the temperature increased. A comparison 
| of the evaporation times for cases (a) 

| and (c) and also for cases (b) and (d) 
shows that t_ decreases. This result 
agrees with our intuition; that is, 

the higher the temperature of the 
surrounding medium, the greater the 


evaporation rate of the droplet. 
Similarly, increasing the relative 
humidity while holding the temperature 
of the medium constant, results in an 
increase in t_. This is seen by compar- 
ing the evaporation times corresponding 
to cases (a) and (b) and to cases (c) 
and (d). It is possible to compare 
other variations of the variables and 
each comparison leads to a result which 
is in accord with our intuition. 


Equations (13) and (14) reveal that, 
for an evaporating stationary droplet, 
the droplet surface area decreases 
linearly with time. There is ample 
experimental evidence to substantiate 
such variation providing the droplets 
are larger than 10um in diameter. See 
for example, Green and Lane. Since the 
diffusion coefficient varies from one 
fluid to another under equal environ- 
mental conditions, the slopes of the 
straight lines defined by equations 
(13) and (14) vary accordingly. Thus, 
in the Maxwell-Langmuir evaporation 
model, the evaporation rate varies 
directly as the diffusion coefficient 
of the evaporating fluid vapor. 


Equations (8) and (9) may be written as 


laM D._ 

Sdt rx Memos, (8a) 
and 

S dt rrt‘P Po (9a) 


These equations show that as the 

radius of the droplet approaches zero, 
the rate of loss of droplet mass per 
unit surface area increases without 
bound. In fact, the loss rate could 
exceed that predicted by statistical 
mechanics. Hence, the Maxwell-Langmuir 
model of evaporation does not properly 
describe the evaporation of droplets 
whose radii become very small. The 
Maxwell-Langmuir model agrees with 
experimental results for droplets whose 
radii are as small as 10um, and there 
is not a large variation from the 
experimental results for droplets whose 
radii are lym or 2pm. 


The next section describes a modifica- 
tion of these equations by Fuchs to 

more properly account for the evaporation 
of very small droplets. 


Fuchs' Modification of the Maxwell- 
Langmuir Model 


The Maxwell-Langmuir model was modified 
by Fuchs who assumed that: 

e The diffusion of vapor molecules 
into the surrounding medium begins at 
the surface of a spherical shell of 


thickness A surrounding the droplet. A 
is a very small distance, approximately 
equal to the mean free path of a molecule 
Geile, 10) 






Spherical 


ke Shell 


Figure 1.--Spherical shell surrounding 
droplet. 
@ The region in the shell is 


nearly vacuous and consequently the 
molecules escape from the droplet and 


‘the region exterior to the shell into 


the shell as if they were escaping into 
a vacuum. 


It is known that the surface tension of 
a very small droplet has a significant 
effect on the rate of escape of molecules 
from the droplet. Surface tension has 

a Similar effect on the rate of adhesion 
of vapor molecules in the surrounding 
medium which strike the droplet surface. 
This effect is accounted for by intro- 
ducing the notion of an evaporation 
constant, a, whose magnitude varies with 
the size of the droplet. Now let C 

be the saturation vapor concentration, 

C_ be the vapor concentration at the 
surface of the spherical shell and € be 
one-fourth of the mean absolute velocity 
of the vapor molecules. Then, from 
statistical mechanics, the mass rate of 
accumulation of molecules from the 
droplet into the shell, Par is 


F_ = 4nr? 
ig abc. 


d ° (18) 


Let F. denote the mass rate of 
condensation of vapor molecules 
originating exterior to the shell 
that strike the droplet, adhere, and 


become part of the droplet mass. Then 
F = AnraEC ° (19) 
s Ss 
In these expressions § is given by 
E | kT (20) 
ptm 


where m is the mass of a vapor molecule 
and k is the gas constant per molecule. 


Thus, the net or observed rate of 
evaporation, Fue is 


= - F 
a a s 


HW 


Zz 2 5 28k 
4nx*ak (C.-C.) (21) 


Now, with the aid of equation (4), the 
mass flux of droplets, F_, escaping by 
diffusion into the medium surrounding 

the shell is 


Pet= 47 le TAD (Cl—-C) (22) 
e co 
where C. is the vapor concentration in 


the surrounding medium. A mass balance 
at the shell surface gives 


aes eee 
n 


© avn Es 


or 
4n(x + A)D(C_-C,) = 4mr7ab(C,-c_) (23) 


Hence, the concentration at the shell 
surface is 


abx*C, + (xr + A) DC. 
oe aer- + (xr + A)D * (24) 


The fuhdamental evaporation equation is 
obtained by equating the time rate of 
change of the droplet mass with the 
rate at which vapor mass is escaping 
into the surrounding medium, Thus, 


ai Powe de (EFA (C_=C,) 
e 


dt 
Or 
2 = 
ath 4n(r + A)Doér (C, C.) 
Og aE ape alk ae 
at abre + (r + A)D » (25) 


From the work of Bradley et al an 
estimate of A is 


Ph em, 4 nm) /m)° (26) 


where m, is\ the mean mass of a gas 
molecule in the surrounding medium, m., 
is the mass of an evaporating molecule 
and X is the mean free path of an 
evaporating molecule. 


Discussion of Fuchs Model 


The mass rate of evaporation per unit 
area is 


Lam _ agD(r + A) (C,-C,,) 
S dt atr-+ (r + A)D : (27) 


As the droplet radius approaches zero, 
the mass rate of evaporation per unit 
area is seen to approach 


B 


aM 

=——s= C.-C 

S at ae ( d = (28) 
which is a finite value. This is in 
contrast to the infinite value approached 
by the mass rate of evaporation per 


unit area given by the Maxwell-Langmuir 
model of droplet evaporation. 


Expressing the mass of a spherical 
droplet in terms of its radius, and 
integrating equation (25) gives an 
equation for the temporal variation of 
the droplet radius. This equation 
permits the calculation of the evapora- 
tion times of droplets of various 
initial radii. Bradley et al. show 
that for small droplets of an organic 
substance, the Maxwell-Langmuir model 
predicts an evaporation time for 10yum 
droplets that is 95.6 percent of the 
evaporation time predicted by the Fuchs 
model. This is not a significant 
difference. In contrast, for droplets 
whose initial radius is 0O.lum, the 
Maxwell-Langmuir model predicts an 
evaporation time that is approxi- 
mately 16.6 percent of that 

predicted by the Fuchs model whose 
results are in close agreement with 
experiment. Thus, the Fuchs model is 
perferred when dealing with droplets of 
very small radii. 


Obtaining the constants appearing in 
the Fuchs model is more difficult than 
obtaining the constants used in the 
Maxwell-Langmuir model. Also, aerial 
spray droplets whose radii are less 
than 5yum are not of sufficient interest 
to be of concern in estimating the 
probability of kill of a target insect 
or pest. However, such droplets are 
certainly of interest in estimating the 
environmental impact of an aerial spray 
Operation. Finally, the analytical 
expressions resulting from the Maxwell- 
Langmuir model are more readily integrated 
into other analyses. Consequently, we 
will not pursue the development of the 
Fuchs model. 


EVAPORATION OF MOVING DROPLETS 
a 


This section of the report will describe 
modifications of the expression for the 
time rate of change of droplet radius to 
account for the motion of the evaporating 
droplet relative to the surrounding 
medium. The presentation will consist 

of a summary of the work of other 
investigators which is deemed most 
appropriate for aerial spray applications. 


The technique most frequently used to 
account for the motion of an evaporating 
droplet relative to the surrounding 
medium is to modify the expression for 
the mass flux of the evaporating mole- 
cules escaping from a droplet assumed 

to be at rest with respect to the 
surrounding medium by a multiplicative 
term of the form 

(29) 


Wee Les £ (rv) 


where v is the velocity of the droplet 
relative to the surrounding medium. 

The function, f(r,v), has the property 
of vanishing when the relative velocity 
is zero, that is, when the droplet is 
at rest with respect to the air. Thus, 
the resulting expression describing the 
evaporation from a moving droplet would 
reduce, in the case of zero relative 
velocity, to the expression describing 
the mass flux of vapor from a droplet 
at rest. Wis frequently termed the 
wiping or ventilation factor. 


Let the "at rest" evaporation rate be 
denoted by ara , then the modified 
equation for the evaporation of a drop- 
let is 


= 


gM _ aM 
dt [2 earn) ] dele (30) 


With the aid of equations (8) and (9) and 
the expression for the time rate of 
change of a spherical mass, the time 

rate of change of the droplet radius is 
given by 


dr 
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(31) 
or 
Dm 


are [2 + £(r,v) | (p-p,,) 


dr 


dt s (32) 


A frequently used ventilation factor is 
that developed by BOS Sisncea (93.8) seu malate 
is 


a5 


W= 1+ a (RE) (338) 


where: 


a is a constant characteristic of 
the evaporating substance 


2vr 
Re = The es Reynold's Number 
v= a = kinematic viscosity of 
s 


the surrounding medium 
H = viscosity of the surrounding 
medium, and 
Rigel density of the surrounding 
medium. 


It is evident that the wiping factor 
increases with the relative droplet 
velocity and the radius. A discussion 
of the effects of Fr&ssling's expression 
on the evaporation rate of droplets is 
given in Green and Lane. They show 
that W increases by about 30 percent 
when the droplet radius is changed from 
5um to 50um. There is considerable 
discussion in the literature concerning 
the numerical value of a. Most of the 
quoted values are such that, 
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Kinzer and Gunn (1951) indicate that 
for droplets whose radii are about 
160um, the value of a is about twice 
that used by Fréssling, while for 
droplets whose radii are less than 
80um, the value of a should be near 
zero. Several investigators, such as 
Lambert and Quince (1969) and Green and 
Lane, set 
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a= 0.276/ (Dp /u)* 


Another form of the Fréssling factor is 


/3 (pny lV? (35). 


Ww=1+ 0.30(sc)* 


where 


Sc = = Schmidt Number 


vI< 


This form is used by Sjenitzer (1963) 
and others as a starting point for 
deriving more complex forms of the 

wiping factor. Pruppacher and Rasmussen 
(1979) present an extensive discussion 

of the determination of the wiping factor 
for large water droplets falling at 
terminal velocity in air. 


Still another technique used to account 
for the ventilation of the droplet by 
the relative motion of the surrounding 
medium consists of altering the diffusion 
constant in accord with the variation 

of the droplet temperature. For example, 
Penski (1978) uses the wiping factor 
expressed in the form given above by 
Sjenitzer, but assumes that the diffusion 
coefficient is a linear function of the 
droplet temperature. Such a technique 
requires the simultaneous calculation 

of the droplet temperature. Several 
expressions for the variation in droplet 
temperature have been presented in the 
literature. A typical expression is 
given by Green and Lane who further 
state that, for droplets whose radii 

are less than 100um, the heat loss by 
radiation can be neglected. They also 
conclude that, for aerosols in normal 
atmospheric temperature, the heat loss 
from the droplet by convection is so 
small that it can be neglected. 


Since the saturation vapor pressure is 
usually a slowly varying function of 
the droplet temperature, it will be 


assumed that the effects of any change 
in droplet temperature on the saturation 
vapor pressure can be ignored. 


None of the above mentioned investiga- 
tions distinguish between a constant 

or a variable-relative velocity. Con- 
sequently, it is assumed that the results 
may be used in conjunction with the 
equations of motion of a falling droplet 
to enable the calculation of the tem- 
poral change in mass of a falling drop- 
Let: 


DISCUSSION 


Determining the evaporation rates of 
droplets comprising an aerial spray 
cloud is a complex task. The carrier 
liquid for a pesticide spray is usually 
oil or water based and the active toxic 
ingredient may be liquid or particulate. 


‘Thus, to properly determine the droplet 


evaporation rate, the rates of evapora- 
tion of droplets composed of multiphase 
liquids and particulate matter should 
be analyzed. 


The spatial distribution of adjacent 
droplets also affects the evaporation 
rate of a single droplet within a spray 
cloud. In this report, for ease of 
analysis, it has been assumed that the 
evaporation rate of a single droplet can 
be estimated by ignoring the influence 
of surrounding droplets. In aerial 
spray applications, this assumption is 
probably reasonable for spray droplets 
far removed from the nozzle. The 
inclusion of the other factors in the 
analysis of droplet evaporation is very 
difficult and most of the existing work 
attempting to include these factors 
results in the construction of complex 
and lengthy computer codes to carry out 
the analysis. 


Other factors affecting the evaporation 
rate are: 

@ Droplet temperature and its 
effect on the density of the droplet. 

e Heat of vaporization, the vapor 
pressure and the surrounding microclimate. 
@ Changes in temperature, atmos- 
pheric conditions, and droplet velocity. 

e Turbulence, high-speed shear flow, 
and fluid viscosity. 
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In addition, the actual field working 
conditions under which spraying is 
conducted are such that spray composi- 
tion, atmospheric conditions, etc., may 
change considerably from one spray 
cycle to the next. 


Finally, other phenomena, such as spray 
release height and swath width, are 
thought to be at least as important, if 
not more important, than droplet evapo- 
ration rate in determining the temporal 
and spatial variation of a spray cloud. 
Thus, it is reasonable to assume that 

a pure liquid evaporation calculation 
may be used as a "first approximation" 
in accounting for droplet evaporation 
in an aerial spray behavior model. 


CONCLUSIONS 


Both the Maxwell-Langmuir model and the 
Fuchs model of the evaporation of a pure 
liquid droplet have been discussed. 
Expressions have been derived which will 
enable the estimation of evaporation 
rates which may be used for comparison 
with experimentally obtained evaporation 
data. The effect of the motion of the 
droplet relative to the surrounding 
medium has been accounted for by 
introducing a "wiping factor." Several 
forms of this factor have been presented 


and the relative merits of each discussed. 


Very small droplets (droplets whose 
diameters are less than 10 micrometers) 
are thought to be ineffective in 
controlling most forest pests. Conse- 
quently, in accounting for the effect 
on insect mortality of the evaporation 
of forest insect sprays, it is appro- 
priate to use the Maxwell-Langmuir 
evaporation rate equation. In contrast, 
the environmental problems associated 
with the drift of aerial sprays are 
known to be caused by very small droplets. 
Thus, when accounting for the effect of 
droplet evaporation on drift, Fuchs' 
evaporation equation should be used. 


Herbicide spraying is usually character- 
ized by droplets whose diameters are 
larger than 500 micrometers. For such 
sprays, the proper accounting of the 
effect of evaporation will necessitate 
the use of a ventilation factor. 


Several experimental programs have as 
their objective the attainment of 
evaporation rate data. The use of such 


‘data in predicting the behavior of an 
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aerial spray will require the use of 
interpolation. The evaporation equations 
derived above should provide the basis 
for an effective interpolation scheme. 
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